Astronomy & Astrophysics manuscript no. main 


©ESO 2013 


March 4,2013 





The evolution of the X-ray emission of HH 2 
Investigating heating and cooling processes 

P. C. Schneider 1 , H. M. Gunther 2 , and J. H. M. M. Schmitt 1 



1 Hamburger Sternwarte, Gojenbergsweg 112, 21029 Hamburg 
e-mail: cschneider@hs . uni-hamburg . de 



2 



Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA, USA 



o 

(N 



in 

6 



> 

in 
o 

On 
cn 

in 
o 

(N 



X 



Received .. / accepted . 



ABSTRACT 



Young stellar objects often drive powerful bipolar outflows which evolve on time scales of a few years. An increasing number of these 
outflows has been detected in X-rays implying the existence of million degree plasma almost co-spatial with the lower temperature 
gas observed in the optical and near-IR. The details of the heating and cooling processes of the X-ray emitting part of these so-called 
Herbig-Haro objects are still ambiguous, e.g., whether the cooling is dominated by expansion, radiation or thermal conduction. 
We present a second epoch Chandra observation of the first X-ray detected Herbig-Haro object (HH 2) and derive the proper-motion of 
the X-ray emitting plasma and its cooling history. We argue that the most likely explanation for the constancy of the X-ray luminosity, 
the alignment with the optical emission and the proper-motion is that the cooling is dominated by radiative losses leading to cooling 
times exceeding a decade. We explain that a strong shock caused by fast material ramming into slower gas in front of it about ten 
years ago can explain the X-ray emission while being compatible with the available multi-wavelength data of HH 2. 

Key words. Herbig-Haro objects - ISM: individual: HH 2 - ISM: jets and outflows - X-rays: ISM - stars: winds, outflows 



1. Introduction 

During their evolution to the main-sequence young stellar ob- 
jects (YSOs) can drive powerful outflows as long as accretion 
proceeds and it is now widely accepted that accretion and ejec- 
tion a re intimately connected dCabrit et al.lll990t lHartigan et alJ 
Il995l) . Protostellar jets can remove angular m omentum from 
the sy stem and thus allow accretion to proceed (Mat t & Pudritzl 
2008). The beautiful manifestations of the outflows are Herbig- 
Haro (HH) objects which are shocked regions wi thin the outflow, 
often also called knots (for an overview see Reip urth & Ballvl 
|2001|) . These HH objects were originally discovered as nebular 
emission lines in star forming regions and this forbidden line 
emission is still a characteristic feature of HH objects. To date, 
HH objects have counter-parts from cm- to X-ray-wavelengths. 
They evolve on time scales of a few years and thus multi-epoch 
observations strongly improve the constraints on the underlying 
physical processes such as the shock heating and the dominant 
, cooling process. 

Many observations of HH objects use optical forbidden 
emission lines (FELs) as tracers for the jet and the tempera- 
ture range exceeding T » 10 4 K is often neglected. However, 
this high energy part of the jet provides important clues on 
the highest velocity part of the jet and thus on the driving en- 
gine which has to be able to provide this high velocity mate- 
rial. The highest temperatures in protostellar jets have been in- 
vestigated with X-ray observations tracing a plasma tempera- 
ture T > 10 6 K. HH 2 cons titutes the first dete ction of X-ray 
emission from an HH object dPravdo et al.l l2.001 ) and following 
this initial discovery, about a dozen protostellar jets have been 
found to emit X-rays. However, most of them are rather dim in 
X-rays and detailed studies of their properties are complicated 
and often ambiguous. Multi-epoch observatio ns exist for two X- 
ray bright sources, the DG Tau jet (HH 158, iGiidel et al.ll2005t 



ISchneider & Schmittf200aiGiidel et al.l200llGudel et al.l201 ll) 
and the L1551 IRS 5 jet (HH 154. IFavata et alJ2002tlBailv et alJ 
120031: iFavata et al.ll2Q06h ISchneider et al.ll201 II) . They show sta- 
tionary X-ray emission close to the driving sources and this 
behaviour greatly differs from what is observed in the opti- 
cal, namely individual knots moving along the jet with proper- 
motion velocities of about 300 km s . However, there have also 
been clai ms of moving X-ray emission regio ns within protostel- 
lar jets dFavata et al. 2006; St elzer et al.ll2009l) . Furthermore, the 
stationary X-ray components are located much closer to their 
drivin g sources than mo st of the optically observed, moving 
knots (Eisloffel &Mundtlll998l) as a resolution of about 0'.' 1 is 
needed to monitor the inner regions of protostellar jets, which 
are also often heavily absorbed. It is therefore a fundamental 
question whether the X-ray emission from HH objects at some 
distance to the driving source shows the same behavior as the 
close-in sources or whether there is a dichotomy with moving 
and stationary X-ray sources in protostellar outflows. 

HH 2 is ideally suited to monitor the time evolution of its X- 
ray emission: It is well isolated so that no contamination from, 
e.g., its driving source is present, it is one of the brightest X- 
ray jets and the first observation dates back to 2000 so that the 
baseline is sufficiently long. 

HH 1/2 

Together with HH 1 , which is located at the opposi te side of the 
deepl y embedded driving source (HH 1/2 VLA 1. iPravdo et al.l 
119851) . HH 2 ranks among the optically brightest HH objects. 
The overview of the HH 1/2 region in Fig. Q] shows that HH 1 
is north-west and HH 2 south-east of the central source which 
is invisible in the optical and near-IR. HH 1/2 VLA 1 is a 
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Table 1. Analyzed observations of HH 2. 

HH 1 



CS star 
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VLA 
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N 



HH 2 
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Fig. 1. Overview of the HH 1/2 region (HST WFPC2 Ha image 
from 1997). The red box indicates the region shown in Fig. |2]and 
the Cohen-Schwartz star is also labeled (CS star). 



class source possessing an in-falling envelope (Reipurth et al. 
119931) . With its loc ation in the Orion s tar forming region at a 
distance of 414pc dMenten etal.ll2007l). multi-ep och HST ob- 
servations (Ba lly et al.ll2002 ; Hartigan et al.j |201 ll) revealed the 
evolution of the individual condensations (knots) within HH 1 
and HH 2 over the course of only a few years. They show that 
the individual knots move with velocities of about 300kms~' 
almost in the plane of the sky. HH 2 appears as highly frag- 
mented with some of the regions showing complex interactio n 
between two knots or reverse facing shocks dHester et al.l ll998). 
A bow-shoc k model for the optica l/UV emission of HH 2 was 
presented by Hartig aiTet al.l (119871) . This model has a shock ve- 
locity of 160 kms~', i.e., some 100 kms~' lower than required 
by the existence of X-ray emission and a post-shock density of 
n e — 2 x 10 3 cm 4 . A comprehens i ve s ummary of the work on 
HH 1/2 can be found in lRaga et all (1201 lh . 

Here, we present a second epoch Chandra observation of 
HH 2 with the aim to constrain the heating and the cooling mech- 
anisms of X-ray emitting HH objects. Inspired by the optical 
data, we assume in the following that the evolution of the X-ray 
emission is rather slow, i.e., that the evolution of the X-ray emis- 
sion of HH 2 is well behaved and not varying erratically so that 
we do not just catch it in two different states by chance. 

2. Observations and data reduction 

The new Chandra exposure uses the same ACIS-S chip as the 
archival observation from 2000, the details of the exposures can 
be found in Tab.Q] We used CIAO 4.3 for the data analysis and 
reprocessed the observations to account for the VFAINT mode, 
removed any pixel randomization and reprojected the second 



Date 


Detector Filter Obs./Prop. ID 


exp. time 




X-ray (Chandra) 






2000-10-08 


ACIS-S 


21 


20 ks 


2010-12-06 


ACIS-S 


12391 


41 ks 




Ho- (HST): 






1997-07-31 


WFPC2 F656N 


6794 


2.0 ks 


2007-08-18 


WFPC2 F656N 


11179 


2.0 ks 



epoch observation such that the brightest sources on the same 
ACIS chip as HH 2 are co-aligned. We estimate relative posi- 
tional uncertainties of about 07 1 based on the remaining offsets 
of 0703, a/08 and 07l4E Spectr a and associated responses were 
extracted using the CIAO tool specextract. 

To compare the X-ray data with the optical, specifically 
Hubble Space Telescope (HST) Ha images of HH 2, we re- 
trieved the corresponding data from the MAST a rchive. The op- 
tical d ata was already discussed in great depth by Hartigan et al.l 
and serves only as a morphological reference for the X- 
ray data here. For the HST data, we relied on the pipeline pro- 
cessing. The alignment of the optical data was performed by 
aligning the Cohen-Schwartz star (V2669 Ori) and checking the 
rotational alignment with the HH objects themselves. 



3. Results 

Our new observation confirms the detection of HH 2 in X-rays; 
the new X-ray image as well as the archival one from 2001 are 
shown in Fig. [2] We find nine photons at the expected position 
of HH 2 using a circle with a radius of 273, while the expected 
number of background photons within this region is only 0.2 
(0.3-1.0 keV) and 0.7 (0.3-3.0 keV), respectively. The exposure 
from 2000 collected 1 1 photons in a similarly shaped region on 
a background of 0.1 (0.3-1.0 keV) and 0.3 (0.3-3.0 keV) pho- 
tons, respectively. The chance probability that any of the of two 
detections is a random background fluctuation is below 10~ 6 . In 
both observations, all photons attributable to HH 2 have energies 
below l.OkeV. 



3.1. The spatial photon distribution 

We calculated the photon centroid using the photons within 
the estimated source region and find that the centroid position 
moved by 171 within the last decade. The position angle of the 
offset is 148° and thus compatible with the optical jet axis. The 
second moments of the spatial photon distributions are 1720 and 
0785 for the 2000 and 2010 exposures, resp.; the spatial profiles 
are shown in Fig. [3] Inclusion of the SE photon visible in Fig. [2] 
for the 2010 observation changes the centroid only by 0/3 and 
the second moment increases to 1709. In order to address the 
significance of these results we used MARX0 simulations and 
two-sided KS-tests. 

For our MARX simulations, we considered (a) a Gaussian 
source flux distribution and (b) a uniformly emitting disk. In both 
cases we used a spatial extent which reproduces the mean second 
moment of the two observations, i.e., cr = 1705 for case (a) and 
a radius of 174 for case (b). Additionally, we investigated a uni- 
formly emitting disk with r = 177 which we regard as an upper 
limit on the size of the emitting region (case c). For the Gaussian 



There is no remaining systematic trend, 
http : // space . mit . edu/ cxc/marx/ 
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Fig. 2. Proper-motion of the X-ray emitting plasma. Red squares 
denote pixel with single photon events and black squares display 
double photon events. Left: X-ray emission (0.3-1.0keV) from 
2000 with Ha contours from 1997 in black. Right: X-ray emis- 
sion from 2010 with Ha contour from 2008. The thick circles of 
radius 2" 3 coincide with the photon centroid, the thin circle in- 
dicates the position of the photons during the 2000 observation. 



model (case a), the surface brightness drops below 10 % of the 
peak value at about 2'.'3; the circles shown in the Fig.|2]have ex- 
actly this size and encircle > 95% of the simulated photons. The 
simulations (a) and (b) show that it is highly unlikely to find the 
measured offset without any proper-motion (< 5%). The 68% 
and 90% errors of the relative offset are 076 and 078, respec- 
tively. The significance of proper-motion decreases to 88 % for 
case (c) and, e.g., the 68 % error increases to 0'.'7. Performing the 
simulations using the centroid calculated from photons within a 
radius of 4" increases the uncertainties by less than 072. 

Comparing the differences of the observed second moments, 
the simulations indicate no significant difference in the underly- 
ing emitting region (< 80% sign.) which is compatible with the 
result of the two-sided KS test giving a probability of 77 % that 
both distributions are not drawn from the same distribution. The 
impact of the off-axis position of HH 2 during the 2000 observa- 
tion affects the result of the KS-test only marginally (the second 
moments change only by < 071). 

In summary, we prefer the conservative interpretation con- 
cerning the evolution of the source shape that no significant ex- 
pansion occurred between the two exposures. On the other hand, 
we regard it as very likely that the majority of the X-ray emitting 
plasma moved in the direction of the optical outflow. 

In Fig.|2]the circles and the lowest contours give the bound- 
ary where the flux reaches 10 % of its maximum value. We note 
that the size scales of the optical and X-ray emission regions ap- 
pear to be similar, with the X-ray emission region being some- 
what larger in both observations. A deconvolution of the optical 
images is not necessary here, as the spatial extent is much larger 
than the point spread function^. 

3.2. X-ray luminosity and energy distribution of the photons 
The observed X-ray count-rate of HH 2 dropped from (3.9 + 
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Fig. 3. Radial distribution of the photons from the centroid po- 
sition with MARX models. The offaxis position of HH 2 affects 
the resulting photon distribution for the extended models only 
marginally (cf. Gauss 2000 and 2010) and the disk models are 
shown only for a single epoch. 



optically thin single temperature plasma (apec) with half so- 
lar abundances. We used cstatistic and binned the spec- 
tra to 22 eV bins. Due to the low number of photons in 
each individual observation, we decided to fix the absorp- 
tion to differe nt values as sugg e sted by the literature, e.g, 
E(B-V)=0.34 (iBrueel et al.l Il98ll iBoehm-Vitense et al.1 Il982t 
iBrugel et al.lll982l) A,/=0.68 dBoehm & Brugellll979T or E(B- 
V)=0.1 1 dHartmann & Ravmondlll984l) ~ 

We fitted the composite and the individual spectra, the results 
for the composite spectra are listed in Tab. The parameters 
for temperature and emission measure (EM) obtained by fitting 
the individual exposures are consistent with the results for the 
composite spectra showing that no significant spectral evolution 
occurred. The fitted temperature is T as 10 6 K while the value 
for the emission measure, and thus the unabsorbed X-ray lu- 
minosity, strongly depends on the assumed absorption (roughly 
EM ~ 10 53 cirT 3 and ~ 10 29 ergs~'). Inspection of the individ- 
ual photons shows that most of the observed emission is consis- 
tent with being O vn emission with a small contribution of O vm. 
This further strengthens the results of the spectral fits as only for 
temperatures around or slightly above 10 6 K such a pattern is 
expected. Although the slope of the effective area changed, the 
mean energy of the photons stays virtually constant with values 
of 543 eV (0.3 keV-3.0keV, 1 1 photons) and 532eV (9 photons). 
This can be explained with the fact that most of the detected pho- 
tons originate from O vn emission (E « 561 - 574 eV), i.e., the 
majority of the photons is emitted within a small energy range. 

In order to estimate which reduction in count-rate is still 
compatible with the data, we calculated the probability for de- 
tecting the observed photon numbers given two different expec- 
tation values, i.e., we calculated 



1.8) x lO^s- 1 to (1.7 + 0.8) x 10- 4 s-' within the last decade. p( n < U,m > 9\m,p 2 ) = CDF(\ x (1 - CDF(%,p 2 )) (1) 



However, the effective area at the very low energies of the pho- 
tons also reduced by about a factor of two during that time 
so that a constant X-ray luminosity is likely. In order to in- 
vestigate the p roperties of the plasma, we used XSPECvl2.5.0 
dArnaudl [T996) to fit the spectra assuming an absorbed (phabs) 



The half energy size of HST is about 07 1 at these wavelengths. 



for different values of p\ and jj 2 using the Poissonian cumula- 
tive distribution function CDF. CDF{x, fi) gives the probability 
to find fewer or equal than x photons given the expectation value 
p. Requiring P > 0.32 (lcr) the ratio p\/p 2 is less than 1.2 while 
requiring P > 0.1 this ratio increases to 1.8. Therefore, the max- 
imum decrease in X-ray luminosity of HH 2 is a factor 1.2 (lcr) 
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Table 2. Results of the specral fits to the X-ray data. 



10 1 



Time (yrs) 



Fig. 4. Evolution of the expected count rate assuming Nh = 
1.2 x 10 21 cuT 2 . The data point shows the measured count rate 
evolution corrected for the changes in the effective area. 



or 1.8 (90% conf.) assuming no spectral evolution as indicated 
above. 



4. Discussion 

Both observations are separated by one decade which is suffi- 
cient to exclude some scenarios for the heating and cooling of 
the X-ray emitting plasma based on the spatial evolution of the 
X-ray emitting region and the evolution of the X-ray luminosity. 

4.1. Proper-motion of the X-ray emission region and its 
spatial morphology 

The new data show proper-motion of the million degree plasma 
in the direction of the optical knot with a nominal velocity of 
Vx-my ~ 210kms~'. The data exclude both, a stationary source 
and velocities exceeding 380 km s -1 , on the 90% confidence 
level. The nominal centroid proper-motion is lower than the ve- 
locities derived for individual regions within HH 2 knot H from 
the optical data (~300kms _1 , iBallv et al. 2002) but considering 
the relatively large error on the X-ray proper- motion, a velocity 
as high as that of the Ha emitting plasma is possible within the 
lcr error. This result for the proper-motion of the X-ray emission 
region of HH 2 represents the most accurate measurement of 
proper-motion of X-ray emitting HH objects to date as there is no 
ambiguity in assigning the initial position of the X-ray emission 
region and there is no contamination from other X-ray sources, 
e.g., from the driving source. The major uncertainties for HH 2 
are the structure of the emission region and the low photon num- 
bers. Longer exposures are required to overcome these issues 
and would also reveal how the X-ray and Ha emission regions 
differ. 

The derived velocity agrees very well with the velocity of 
« 200km s" 1 of the outer X-ra y emitting jet of DG Tau (T x 
3.4 x 10 6 K, iGiidel et alJl201lh . On the other hand, the veloc- 
ity of HH 2 is con siderably smaller th an the X-ray emitting knot 
in the Z CMa jet dStelzer et al.ll2009l) . However, the initial posi- 
tion of the latter knot is rather uncertain and the estimated X-ray 
proper-motion exceeds the optical estimate in that case (contrary 
to HH 2 and DG Tau). We do not include the jet of L 155 1 IRS 5 
here because the majority of the X-ray emission from this jet 
is in fact stationary and no reli able proper-motion value can 
be derived dSchneider et al.l201 ll) . Considering the temperatures 
of the X-ray plasmas, it appears that higher post-shock proper- 
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motion leads to higher plasma temperatures. A similar trend has 
been found for individual knots in HH 2 between the proper- 
motion velocity and the line width which is an in dicator of the 
actual shock velocity by Eisloffel & Mundtl(ll998l) . However, we 
refrain from claiming a trend because the number of sources is 
low and because the post-shock velocity is not necessarily re- 
lated to the actual shock- velocity. 



4.2. Heating and cooling 

In principle, there are two possibilities for the constant X-ray 
luminosity: Either the cooling time is rather long or the X-ray 
emitting plasma is constantly replenished. By investigating the 
importance of adiabatic, radiative and conductive cooling we ar- 
gue that slow cooling is the most likely scenario and discuss the 
resulting implications for the heating process. 



4.2.1. Cooling 

The X-ray emission region of HH 2 is incompatible with a point 
source in both exposures, a s already shown from the first epoch 
data bv lPravdo et al.l(l2001l) . Furthermore, we do not find any ev- 
idence for an increasing (projected) size of the emission region. 
Assuming a sphere with the same radius as the circle in Fig. |2] 
and a volume filling factor of unity, we derive electron densi- 
ties ofn e ~ 10 2 - 10 3 cm depending on the exact value for the 
emission measure. For smaller emission regions, i.e., if the emis- 
sion is composed of individual clumps, the electron density will 
increase. With these densities the X-ray emitting material can be 
in thermal pressure equilibrium with the other temperature com- 
ponents and the absence of significant source variability makes 
cooling by expansion unlikely. 

Using the CHIANTI atomic database jDere et alJ 1 19971 

2009), we estimate that the observed count rate would drop 
by a factor of 1.4 within 10 years by radiative cooling assum- 
ing n e — 5 x 10 2 cm- 3 or by 20% within 30 years assuming 
n e - 100 cm- 3 (ignoring any further drop in the effective area, 
see Fig. @), which is perfectly compatible with the data in the 
absence of reheating. The variation in the temperature is even 
smaller: For n e = 10 3 cirr 3 the temperature changes by only 
10 % within a decade so that the appearance of a constant plasma 
temperature is fully compatible with radiative cooling. 

Using the approximation given in Orlando et al ] (120051) for 
the conductive cooling time of a million degree plasma embed- 
ded in a gas of much lower temperature, we find that the con- 
ductive cooling time exceeds 10 4 years in the case of HH 2 and 
is thus negligible compared to radiative cooling. 

In summary, the dominant cooling process is most likely ra- 
diation as the cooling time of thermal conduction is too long, 
while the data does not indicate cooling by expansion. 
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4.2.2. Heating 

Currentl y, knot H is in a phase of being overtaken by the faster 
knot M dEisloffel et al.lll994t iHartiganet ail 1201 ll) , so that the 
X-ray emission might not be intrinsic to knot H but rather to 
the interaction of both knots. The me asured velocities fo r the 
two knots were 220 and 490 km s _1 (Eisloffel et all 1 19941) re- 
sulting in a shock velocity of v S h = 270 km s -1 . The derived 
temperature of the X-ray emitting plas ma of about T « 10 6 K 
requires shock velocities of 260km s -1 dRaga et al] |2002). i.e., 
the expected shock velocity is precisely the value needed for 
the observed X-ray emitting plasma. The expected post-shock 
proper-motion velocity would be 290kms~' which is compati- 
ble with the X-ray and Ho- data. Even the nominal differences 
between the proper-motion of X-ray and Ha emission would 
not lead to detectable offsets during this time span. This picture 
can also explain the non-detection of O vi emission (peak for- 
matio n temperature 3.2 x 1 5 K) in the HUT spectrum taken in 
1995 (Raymon d et alj 19971) by assuming that knot M did not en- 
counter knot H by that time and the existence of X-ray emission 
in 2001 would be a fortunate event in the history of HH 2. 

Concerning the properties of the outflow, this scenario re- 
quires a jet with vastly different velocity components which 
move approximately into the same direction. As the flow-times 
of the two ejecta to the current position differ, both knots result 
from different outflow epochs. However, it is not clear whether 
knot H is regularly fed by a high-velocity, pulsed outflow or 
whether these events are rather stochastic and less frequent. 
Since knot M was observable before encountering knot H, it also 
experienced shocks before. These shocks, however, have to be of 
much smaller shock velocity, first, because knot M preserved its 
high velocity and, second, because no O vi emission was de- 
tecte d from the correspon ding region in the 1995 HUT observa- 
tion dRavmond et al.ll 19971) . 

The second possibility to explain the constant X-ray lumi- 
nosity and temperature is to assume that the million degree 
plasma is constantly refueled by a number of smaller shocks ap- 
proximately uniformly distributed over the time. Although the 
interaction between knot M and knot H p ersists for some time , 
the dramatic changes of the Ha emission (Hartigan e t al.ll20l"ih 
argue against a quasi-static shock scenario for the X-ray emis- 
sion. Furthermore, the data does not indicate higher photon en- 
ergies closer to the centroid position as one would expect if the 
X-ray emitting plasma is constantly replenished at the position 
of the Ha producing shock followed by expansion and cooling 
of the post-shock plasma. 

The available X-ray data of HH 2 therefore suggest that a 
single strong shock about a decade ago caused the heating to 
X-ray emitting temperatures since constant re-heating is not re- 
quired, i.e., it is very likely that individual strong shocks heat a 
fraction of the material to X-ray emitting temperatures. 

5. Summary and conclusion 

As protostellar jets and consequently also HH objects are dy- 
namic objects evolving on time scales of years, our second epoch 
observation of HH 2 allows to further constrain their heating and 
cooling processes. The new data show that the X-ray emission of 
HH 2 moves into the same direction as the optical knots. The X- 
ray proper-motion velocity is lower but still compatible to that 
of the Ha emission. The constant X-ray luminosity, the absence 
of strong expansion and the unimportance of conductive cooling 
imply that the cooling of HH 2 is dominated by radiation with a 
rather long radiative cooling time of a few decades. 



With these facts, the most likely scenario is that the X-ray 
emitting part of the knot results from a single, strong shock when 
a bullet of high speed material rams into much slower material 
in front of it. For HH 2, these two components are probably the 
high velocity knot M and the slower knot H since their veloc- 
ity difference is exactly the shock velocity needed for the ob- 
served X-ray emission. There is no evidence for a continuous 
or significant replenishment of hot, X-ray emitting plasma but 
a significant fraction of the optically observed gas might result 
from additional slower shocks. This naturally explains that both 
components are (almost) co-spatial and move with comparable 
velocities into the same direction. 

In order to measure any separation between the X-ray and 
Ha emitting part, one needs longer exposures resolving the 
structure of the X-ray emission. Future observations will also 
reveal the true X-ray cooling time which we could not address 
with our second epoch observation due to the lack of significant 
cooling within the last decade. 
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